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Current-phase relation in a Josephson junction formed by putting two s-wave superconductors 
on the same edge of a two dimensional topological insulator is investigated. We consider the case 
that the junction length is finite and magnetic impurity exists. The similarity and difference with 
conventional Josephson junction is discussed. The current is calculated in the semiconductor picture. 
Both the 2-7T- and 47r-period current- phase relations (Iiir{4')> ^tWO) are studied. There is a sharp 
jump at cj> — n and <j) = 2tt for and 1^ respectively in the clean junction. For I-z-k, the sharp 
jump is robust against impurity strength and distribution. However for I^, the impurity makes 
the jump at <j> = 2ir smooth. The critical (maximum) current of Ii^ is given and we find it will be 
increased by asymmetrical distribution of impurity. 
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I. INTRODUCTION 

Recently the topological insulator (TI) has excited 
great interest in the condensed-matter community.^] 
The unique feature of TI is the existence of edge states 
(or surface states) which is protected by time reversal 
symmetry. The edge state of a two dimensional (2D) TI 
can be considered approximately as a ID mental. But 
since spin and momentum direction of carriers is locked 
together owing to strong spin-orbit coupling, it's only half 
of the ordinary electron gas. This helical property is ro- 
bust against nonmagnetic impurity due to its topological 
origin. If the edge state is in contact with a superconduc- 
tor, a topological superconducting edge stat e will form in 
the interface because of proximity effect And it can be 
viewed as a ID topological superconductor (TS). There- 
fore it's able to construct a Josephson junction on one 
edge of the 2D TI. 

Experimentally, the edge state in HgTe/CdTe quan- 
tum wellsP in InAs/GaSb quantum wellspl and surface 
state in -B^S^ systems^ have been observed. The su- 
perconducting proximity effect and Andreev reflection in 
InAs/GaSb quantum wells and Bi 2 Se3 systems coupling 
to superconducting electrode have been demonstrated!^ 

The conventional superconductor-normal metal- 
superconductor (SNS) junction has been investigated 
in detail in the last three decadesP^HHI Since the 
superconductor-TI-superconductor (STiS) junction is 
only half of the SNS junction, the corresponding Andreev 
bound state^ and current-phase relation are similar for 
the clean junction if we suppose quasiparticles distribute 
thermodynamically (27r-period current case) However 
for STiS junction, a 4-7r-period current-phase (I^ (</>)) 
relation (fractional Josephson effect) may arise if the 
thermodynamical distribution is partially destroyed 
while super cond ucting phase difference is changed 
adiabaticallyPE^ The effect of nonmagnetic impurity 
and magnetic impurity is identical for the SNS junction 
due to spin degeneracy. However for the STiS junction 
only magnetic impurity can lead to a backscattering 



owing to time reversal symmetry. In dirty junctions 
magnetic impurity contributes another significant differ- 
ence, the extra n phase shift for hole reflection!^ As a 
result even the 27r-period current (hiri'P)) an d Andreev 
bound states of STiS junction would be quite different 
from those of the SNS junction. 

However in earlier workpEfiHUl only short STiS junc- 
tion (junction length L far less than the superconductor 
coherent length £o) is studied. And it's only very recently 
we notice that the work by Beenakker et aZp2 discusses 
the finite length clean junction. To the best knowledge 
of us, a study of the finite length STiS junction affected 
by magnetic impurity is still missing. That is the gap we 
want to fill here. 

In this article both the 2-7r-period and 4-7r-period 
current-phase relation is calculated. There is a sharp 
jump at <f> = 7t and 4> = 2tt for 7 2 ir and Ia-k respectively 
in the clean junction. For l^, the sharp jump at <j> = it is 
robust against impurity strength and distribution. How- 
ever for 7 47r , the impurity makes the jump at </> = 2tt 
smooth. The critical current and shape of current-phase 
characteristics are greatly influenced by junction length. 

The rest of the paper is organized as follows. In Sec. 
|ri| we describe the model and give the analytical results. 
In Sec. |III[ the numerical results and analysis are given. 
In Sec. |IV| we give a brief conclusion. In Appendix A, we 
give the reason of the similarity between STiS junction 
and conventional SNS junction. In Appendix B and C, 
we derive the current operator and give the detail of the 
calculation. 



II. MODEL AND ANALYTICAL RESULTS 

Two s-wave superconductors are in intimate contact 
with one edge of 2D TI. Because of the proximity ef- 
fect, a ID TS forms in the interface. Then we have a 
STiS Josephson junction on one edge of the 2D TlP^I 
The effective Hamiltonian of the edge state is given as 
Hq = vpa^Px, in which p x = —ihd x , 04,2,3 are Pauli ma- 
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trices acting in the spin space and vf is the velocity of 
the edge statesP Proximity effect contributes a paring 
term, then the Hamiltonian of the ID TS is given asP 



H 



dxi]y{Ho — fi)ip - 



+ A> ; Vt (!) 



in which ip = (ijj-f, V^) T \ ip-f (^4-) annihilates a right (left)- 
moving electron. A = A e J * is the paring potential, 
A = |A| and <j)' is the phase of the superconductor. 
In Nambu representation ^> = (ifo,ipi, ipp — ip\) T , with 
ihd t ^ = Hbclg^ we d erive the Bogoliubov-de Gennes 
(BdG) HamiltoniarPMl 



H 



BdG 



vfPx&3T 3 - fir 3 + A o [cos(0')ti - sin(<f)')T 2 ],(2) 



where fi is the chemical potential and Ti.2,3 are Pauli 
matrices mixing the ip and ijp blocks of Parti- 
cle hole symmetry is expressed as {HbcLG, 3} = 0, 
in which EE = (J2T2K and K is the complex conjuga- 
tion operator. As a result these states are not inde- 
pendent. For an infinite TS the dispersion relation is 
e = ±y / h 2 Vp(k ± hp) 2 + |A| 2 , in which /i = hvpkp. 
And we neglect the self-consistency condition of AP3 
For the junction considered here, A = Aoe l ^ 1 0(— x) + 
A e l ^ 2 8(x — L) where L is the length of the junction. 

We include a region with magnetic impurity by adding 
a scattering term Ma\6{x — L\)0{L2 — x) in HsciG- The 
magnetic impurity can change the direction of particles, 
which can be described by the scattering matrix for elec- 
trons and holes 



S e = 
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-r*t/t* 



and Sh — 
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We denote the reflection coefficient R = \r\ 2 and transi- 
tion coefficient T = \t\ 2 . For simplicity, we have assumed 
that R is a constant independent of energy and the length 
of the impurity region. Under this assumption the effect 
of the length of the impurity region is equivalent to re- 
placing the junction length L with an effective length 
L' = L — (i 2 — and in the following we abbrevi- 
ate L' to L. Comparing with SNS junction,^ there is 
an extra it phase shift for hole reflection, and that's the 
origin of difference between STiS and SNS junction in 
Andreev bound states and I2K (see Appendix A for an 
explanation) . 

Incident particles with energy e will be reflected at the 
superconductor-normal interfaced For SNS junction, it 
can occur both the Andreev and normal reflections at the 
interface. But for STiS junction, only the qu antum An- 
dreev reflection occurs at the interface! 8 1 31 1 32 1 If |e| < A , 
incident particles will be reflected completely, therefore 
Andreev bound states will formPSl Solve the BdG equa- 
tion, then we obtain the energy level equation of Andreev 
bound states. For clean junction 
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where £0 = hvp/{2Ao) is the superconducting coher- 
ent length, (f> — (f>2 — 4>i is the phase difference and 
n = 0, ±1,±2, ... . The second term on the left side of 
Eq.Q is equal to (k e — kh)L, where k e (kh) is the wave 
vector of the right-moving electron (left-moving hole) 
with energy e. Then we can interpret Eq.Q in terms of 
Bohr-Sommerfeld quantization of the periodic electron- 
hole orbits in the TI region! 3 ^ In the presence of impurity 
Andreev bound state is given as 



, e . e L 
— 2arccos( —— ) + — - — = a 

in which the phase difference is changed to a, 



cos(a) — Tcos(<p) — Rcos{ 



e L-2L X 



(5) 
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which is different from that of the SNS junction.^ 

The Josephson current /(</>) induced by the supercon- 
ducting phase contains two parts, the discrete current 
Id{<P) and the continuous current I c (<f>) carried by quasi- 
particles occupying Andreev bound states and continuous 
energy spectrum respectively. To compute the current, 
we suppose the system is nearly in thermodynamic equi- 
librium. Because the current is constant, we can solve 
the wave function and then obtain the average value of 
current operator in the TI region. The current due to the 
scattering state (the eigenstate of junction Hamiltonian) 
ip = (u(x), u'(x), v(x), v'(x)) T with eigenvalue e is 



J = €Vp\ 
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where e is the electron charge and /(e) is the Fermi distri- 
bution function. The last two terms describe the current 
carried by the "vacuum" (spin-down band and spin-up 
band filled by electrons) on which we can create quasi- 
particles occupying the ground state of HsdG to obtain 
the superconducting ground stateP^ There is an alterna- 
tive statistical method by which current is the derivative 
of free energy. In this article we use the wave function 
method to calculate the continuous current and the quan- 
tum statistical method for the discrete current. In ap- 
pendix B and C we give the calculation detail and prove 
results according to both methods are equivalent for the 
discrete current. 

The discrete current can be written as Id{4>) — 
J2 n In{4>) i ( e n) j where I n (<j>) is the current carried by the 
quasiparticle occupying Andreev bound state with eigen- 
value e n . According to the quantum statistical method, 
the effective current due to Andreev bound state with 
eigenvalue e„ is I n {4>) — f^J (derived in Appendix C). 
For dirty junction, 



In{4>) = 
7 = 1 
x sin( 



1 evp Tsin{4>) 1 
2L + 2£(e„) sin(a) 7 

h evp R L — 2Li 

2eA L + 2f(e„) sin(a) £ 
e n (L-2Lt). 



A £o 



(8) 



3 



For clean junction, 
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(9) 



is the energy dependent coherent 



where £(e) = 
length. 

For a short junction (L << £ ), it's enough to consider 
discrete current only, because the continuous current is 
of the order of L/£o- However for a long junction the 
continuous current can not be neglected. To calculate 
I c (<j)), we first construct the scattering state for an inci- 
dent particle having energy e, and then apply the current 
formula given byEq.Q. And we take the semiconductor 
picture (both the positive and negative solutions of BdG 
equation are used) . The detail of constructing scattering 
states and computing current is similar to Ref ill I and 
some detail is given in Appendix B. Results are given 
below. For clean junction, 

r—Ao rOO 
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For dirty junction, 
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in which 
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FIG. 1: (Color online) Andreev bound states with several 
junction length and transition coefficients. Red line T=l, 
blue line T=0.5. Left: L = 0. Middle: L = 8£ , LI = L/2. 
Right: L — 8£o, LI = 0.2L. 



of a;. The symmetrical impurity (Li = L/2) opens a gap 
at (f> = 2rvK (n is the integer) as shown in the middle panel 
of Fig.l. For the asymmetrical impurity (L\ ^ L/2), it 
can open the gap at both <f> = 2n7r and <fi — (2n + 1)tt 
(see the right panel of Fig.l). But the crossing point at 
= 7r, e = remains for arbitrary length and can not 
be broken by impurity scattering which is different from 
the conventional SNS junction!^ That specific crossing 
point is protected by the fermion parity conversion!^ 
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FIG. 2: (Color online) Current-phase relation for several junc- 
tion length with certain transition coefficient. T=l (0.5) for 
solid (dotted) line. For dotted lines the impurity distributes 
symmetrically. Current in unit of 2 (l+2^ ) ■ Length in unit of 
£o- Temperature is zero. 



III. NUMERICAL RESULTS AND ANALYSIS 



For the short junction (L << £o), only a pair of An- 
dreev bound states contributes to current, and analyti- 
cal result is available!^ Andreev bound states are given 
as e = ±e , e = >/TAqcos(4>/2). The corresponding 
current is / = I tanh( 2k e ° T ), Iq{4>) = ^ y/TA sin(<fr/2) 
where ks and Tb are Bortzmann constant and tempera- 
ture respectively. However for the finite length junction, 
we mainly give numerical analysis. 

FigJT] shows the effect of junction length and impu- 
rity reflection on Andreev bound states. The length of 
the junction will increase the number of bound states 
consistent with the usual ID quantum wells. The num- 
ber can be given approximately as 2 Int(L / (£qw)) + 2 or 
2 Int(L / fair)) +4, where Int(x) means the integer part 



The zero temperature current-phase characteristics for 
different junction length and impurity strength and dis- 
tribution are shown in Fig(2]and Figj3j With the length 
increasing, the curve changes from sinusoidal to saw- 
tooth. Impurity reflection is mainly to decrease the crit- 
ical current. There is a robust sharp jump at <j) = w. 
Since the continous current is zero while <f> = 7r, the jump 
is rooted in the crossing point of Andreev bound state 
at <j) = 7r, e = 0. It will not be destroyed by impurity 
reflection because the impurity can not open a gap at 
= tt, e = 0, which is different from the case of conven- 
tional SNS junction. 

The critical (maximum) current I c .2tt is reached when 
<f> = 7T, with / Ci 27r = I<i{tt) due to I c (^) — 0. For the dirty 
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For the clean junction = 



2£i) 2 

For symmct- 
In 



2 L+2f 

rical impurity distribution, I c ,2-k\t = VTI c ^tt\t=i- 
this case, the impurity reflection monotonously decreases 
the critical current. The asymmetrical impurity distribu- 
tion will enhance the current shown in the inset of Fig(3] 
That's different from the conventional SNS case where 
the critical current will decrease when impurity leaves 
the center. For the long junction with extremely asym- 
metrical impurity distribution (L » £ , L » Li), we 
have I c .2tt ~ I -L+2g ^ or ^ n0 ^ * 00 sma U, which nearly 
reaches the result of clean junction. 
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FIG. 3: (Color online) Current-phase relation for several LI. 



Current in unit of 



LI in unit of L. T = 0.5, L = 



2(£+2£ )' 

8£o- Inset shows the dependence of critical current on LI. 
Temperature is zero. 



jump located at <f> = 2n smoother for 1^. The reason 
is that for clean junction the energy crossing of Andreev 
bound state at <f> — 2ir has a none-zero slope. While for 
the dirty junction the slope is zero (see Fig[T]). Here we 
denote the maximum of (h-n) as 7 c ,47r = glc.2-n- 9 
increases with length increasing. We have g — 1 for the 
junction with length L = 0. For the long clean junction 
(L >> £o)> g = 2p3 That's apparent if we notice that 
the energy level located deeply in the paring potential 
well is nearly linear for the long clean junction. Impurity 
reflection will make the factor decrease. For long junc- 
tion case, vary the reflection coefficient from to 1, g 
changes from 2 to 1. For a short junction (L = 0), g is 
independent of reflection and we have g = 1. 
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FIG. 4: (Color online) Current-phase relation to show the 47r 
period. Dotted (solid) line for 1^ (h^)- Left: T=l. Right: 
L=8£o, L\ — 0.5L. Current in unit of 2 (l+2£ ) ■ L en gt.h in 
unit of £o- Temperature is zero. 



In the previous discussion, we suppose that there is 
some mechanism to make quasiparticles distribute nearly 
thermodynamically. Now we consider the case that the 
necessary mechanism is absent for the two eigenstates 
ip±{<j)) with energy e±(</>) nearest to zero shown in FigJT] 
The two states are connected by electron-hole symmetry, 
ip + = and e_ = — e+. 

The original state remains while phase difference is 
changed adiabatically. Starting from ground state while 
4> = 0, for tf> < 2ir state e_ is occupied. The current due 
to a pair of Andreev bound states is 

SW* 1 -* ' (14) 

and the distribution is /(e-) = 1 independent of energy, 
then we have I e = | for < <f> < 2tt. While <j> = 2n, 
the system is in excited state. And it can not decay to 
ground state because of fermion parity conversionl^SlFor 
2tt < <f> < 4ir, the state e + is occupied. While <p = An, 
the system reaches the original state we start withpESl 
Therefore I e is 4-7T periodic. The net current will be 4-7T 
periodic since I e contributes significantly to current. 

The current-phase curve is shown in Fig|4j There is a 
sharp jump at <p = 2tt for 1^ in finite length clean junc- 
tion. For 7 2 7rj the jump at 4> — ir is robust against impu- 
rity reflection. However impurity reflection will make the 



IV. CONCLUSION 

In summary, the current-phase relation of a finite 
length STiS junction with magnetic impurity is investi- 
gated. We consider both the 2tt- and 47r-period case. 
With the length increasing, the current-phase curve 
evolves form a sinusoidal shape into sawtooth shape. 
There is a sharp jump at (f> = 7T and (f> = 2ir for and 
Iin respectively in the clean junction. For I2 W , the sharp 
jump at <p = 7T is robust against impurity strength and 
distribution. However for I^, the impurity makes the 
jump at <j) = 2tt smooth. The critical current is greatly 
influenced by junction length and impurity. 
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Appendix A 

This appendix explains the origin of the similarity and 
difference between the STiS junction and conventional 
SNS junction. 

For the SNS junction, Nambu basis can be selected in 
two equivalent form, tp + = (ip^,ip^) T or tp_ — (ip^,ip^) T , 
due to spin degeneracy. With ihd t ip± = Hb<ig'4>±, we can 
derive H BdG = (pl/2m- fi)r 3 + A [cos(<//)tl -sm(^')r 2 ], 
p x = —ihd x and m is the eff ective mass of electron. Take 
the Andreev approximation 30 ! 36 ! anc j denote the eigen- 
vector as tp — ^ e l ( crk F+ Sk ) x ^ a = ± for incident particles 
with wave vector near ikp, x is a vector independent of 
x. Then we arrive at the Andreev equation^ 



cvfPx — hvpkF 



-avpPx + hvpkF 



(p = etp (15) 



If we reset Nambu basis as W = (ipf, ip±, ipp V-'|) T an d 
take a = 1(— 1) for tp + (ip_), we will find the correspond- 
ing BdG Hamiltonian is identical to the BdG Hamilto- 
nian of STiS junction, but the Nambu basises are con- 



nected with a unitary transformation P = 



03 



, m 



which (To is a 2 x 2 unit matrix. P matrix leads that for 
dirty STiS junction there will be an extra 7r phase shift 
for hole reflection as is shown in Eq. ([3| . The other choice 
to take a = 1(— 1) for ip_(i/j+) corresponds to the same 
junction formed on the other side of the 2D TI. 



Appendix B 

This appendix is to derive the current formula Eq.Q 
and give some detail of calculating the current. 
The system is given as 



H = J dx^(H - h^j + A^I + A*^ (16) 

With ihd t ^ = HsdO^, the BdG Hamiltonian is yielded 
H -n A 



H 



BdG 



A* -T{H - n)T~ x 



(17) 



with time reversal operator T = —ia 2 K. In fact Eq.(17) 



is appropriate for arbitrary Ho but with the correspond- 
ing time reversal operator for different systems. The BdG 
equation can be written as 



H 



BdG 



(18) 



where <p iiV = (u i ^(x),u' lu (x),v i ^(x),v' iv (x)) T is the 
eigenvector and is the eigenvalue. Because of 
electron-hole symmetry {PtBdGi^} — 0, l Etpi jU {x) is also 
an eigenvector with eigenvalue —£%, v . v and i denote en- 
ergy and the extra degeneracy respectively. For continu- 
ous spectrum i = 1,2,3,4, tp 3 v = Ep 2 ,u, <fi,v = — S<pi„. 



ipi :l/ is the scattering state constructed from the inci- 
dent state (eigenstate of ID TS) (pi >L , shown in Fig{5] 
However for Andreev bound states we only have i = ±, 
(fi- M = S<^ +! „. For simplicity we denote ip~^ = <pi, v , 
p+, v = <p>A V and ip 2 ,u = <P3,v = 0. 




wave vector 



FIG. 5: (Color online) Continuous energy spectrum for infi- 
nite ID TS. Red and green lines correspond to electronlike 
and holelike eigenstates of ID TS respectively. (f>i V is the the 
incident state from which one can construct the scattering 
state tpiv. 

To diagonalize the Hamiltonian we first rewrite it 
as H = h J dx^ PtBdG 1 ^ + constant. The Bogoliubov 
transformation is given as f = Tli U S l j^ Ul in which 

S v = (<Pl,v,<P2,v,<p3,v,<P4,u), lv = (71,^72,^73,^74,^, 

74. v = — 7J v , 73, „ = 72 „. Jiu annihilates a quasiparticle 

in eigenstate ip iv . Then we have H = ~ £\ w e i^l] M lt.y + 
constant. 

The current density operator can be derived with the 
current density conversion equation, dtp{x)+d x J{x) = 0, 
in which electron density operator p(x) — e(ip^(x)tp^(x) + 
ipl(x)ipi(x)). In the TI region, it can be derived as 

J(x) — evF(ip\{x)^{x) — ip\(x)ipi(x)) . In the TS region, 
the paring potential will contribute an additional term 
—d x J s = 2e(Aiplip\—H.C.)/ih, which describes exchang- 
ing Cooper pairs between quasiparticles and condensate. 
However this term vanishes for energy larger than paring 
potential, thus it makes no contribution to the continuous 
current. But it will make the discrete current transforms 
into supercurrent carried by the condensate gradually in 
the superconducting region.^ 

Take the ensemble average J(x) — <J(x)>, with Bo- 
goliubov transformation and <%^f%,v> = f In the 
TI region we find J = J2 V Ji,v + Ji2,u, 

Ji,v= ev F [\ui^\ 2 f{ei^) ~ K„| 2 (l - /(£»,»/))] 

[|<J 2 /(^)-|<J 2 (l-/(^))] (19) 

which is just Eq.Q we want to derive. The extra 
current owing to paring potential is < —d x J s > = 

1 E MAE i=1 , 2 ,K/<,,/(v) : ^>y(£i,,)]}. . 

Now we prove that the contributions from electronlikc 
and holclikc injected states are equal. J ev (Jh,v) is the 
current due to electronlike (holelike) state tp e>v [}Ph,u) 
with eigenvalue e Ci „ (eh, v ) where e = {1,2}, h = {3,4}. 



G 



J e /h,u is given by Eq.(T9j). Since ip h<v = Ecp etU , e h>v = 
—e e ,v and 1 — /(e) = /(— e), we can obtain J 6il/ = Jh.v 

For continuous spectrum, the eigenstate with a cer- 
tain energy is 4-fold degenerate. The continuous cur- 
rent can be written as, I c = J de v N{e v )J(e v )/2, 
J (tv) = E t J i{tiv) = ev F[J2i(\ u i,v\ 2 + \ v iA 2 - KJ 2 - 



)f{ti,v) ~ Ei I 



^l 2 + E>LI 



! ], N(e v ) is the den- 



sity of states of TS. Solve the eigenvectors and we find 
the last two terms cancel with each other. Then we have 
J c = ~\ I d^uN(€ v )[Ji(e v ) + J' 2 (e v )]tanh( 2k c » Ta ), where 



J[{e v ) = ev F (\u^\ 2 + K,| 2 - |u{ J 2 - |<J 2 ), which is 



2k B T B ■ 
I |2\ 



the equation we use to derive Eqs. ( 10 )-( 12 1 



Appendix C 

This appendix is to prove the discrete current ob- 
tained by wave function method and quantum sta- 
tistical method is identical if the states are occupied 
thcrmodynamically!^ 

In this appendix we take the Nambu basis given as 
vj/' = (xpf, V>J, ip±, — "0|) T for simplicity. The correspond- 
ing BdG Hamiltonian is 



H' 



BdG 



( (VfPx ~ M) r 3 + A 

V M(x)t 



M(x)t 
(-v F p x - h)t 3 + A 



(20) 



where M(x) 
unit matrix and A 



M6[x - Ly)Q{L 2 - x), r is a 2 x 2 

sgn(x) = x/\x\. 

A pair of Andreev bound states connected by electron- 
hole transformation is given as ip'± with energy e±, 
= Eip'_, <p'± = (u ± (x),v ± (x),u' ± (x),v' ± (x)) T and 
e_ = — e+. The corresponding current is J = evp(\u_ | 2 + 



\v-\ 2 



'_| 2 -K| 2 )(/(e-)-/( £+ ))/2 + e UF (| M _h 



,/ 12 



)/2. With the solved eigenvectors, 



we find the second term on the right side vanishes. Then 
the current is derived as 



(21) 



,' 1 2 



,/ 12 



)/2 



J = J( e _)/( e _)+7(e+)/(e + ), 

in which I(e±) = evp(\u±\ 2 + \v±\ 2 
and I(e + ) = — 7(e_). I{e±) can be seen as the effective 
current carried by eigenstate tp±. 

Rewrite the current as J = 7(e_)(f(e_) — /(e + )), 
and then act the operator p x = —ihd x on both 
sides. With a straightforward calculation, we have 

^PxJ/e = <p1 ( [A ' T3] [AjT3] )^'- ' (/(e-) - 

With [A, r 3 ] = 4isgn(x)^, we derive ^sgn(x)d x J = 

V- dH a B / a y/_(/(e_) - f(e+)). Integrate among the the 
whole region. Since the current in the TI region is 
constant and it decays to zero gradually in the super- 
conductor, the left side gives HJ(x = 0)/e. With the 
help of Feynman-Hellmann theorem the right side gives 



<), 



-(/(e_) - /(e + )). Then we obtain 



J 



e de- 
hlkb 



f(e-) 



e de + 
h~d0 



(22) 



Comparing with Eq. (J^lj) we have I(e±) = |^^- So far 
we have proved the two methods are equivalent for the 
discrete current. 

It's of importance to point out that the quantum sta- 
tistical method has taken both particle energy levels and 
hole energy levels into consideration. For a clean STiS 
junction with length L — 0, one can take the Nambu basis 
to be -0 = (4>^ , ipj) T and the corresponding BdG Hamil- 
tonian is 2 x 2. In this case there is only one energy level 
with energy e_(0) contributes to current. However as we 
have discussed a wrong result, J = | ^-/(e_), will be 
derived if we use the quantum statistical method . 
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